Rotating helical turbulence. Part I. Global evolution and spectral behavior 
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We present results from two 1536 3 direct numerical simulations of rotating turbulence where both 
energy and helicity are injected into the flow by an external forcing. The dual cascade of energy 
and helicity towards smaller scales observed in isotropic and homogeneous turbulence is broken in 
the presence of rotation, with the development of an inverse cascade of energy now coexisting with 
direct cascades of energy and helicity. In the direct cascade range, the flux of helicity dominates over 
that of energy at low Rossby number. These cascades have several consequences for the statistics 
of the flow. The evolution of global quantities and of the energy and helicity spectra is studied, 
and comparisons with simulations at different Reynolds and Rossby numbers at lower resolution are 
done to identify scaling laws. 



I. INTRODUCTION 

Helicity, the alignment of velocity and vorticity, is a 
measure of the number of links in vorticity field lines, and 
an indication of lack of mirror-symmetry in a flow 
Isotropic and homogeneous turbulence with and with- 
out helicity has been thoroughly studied in the literature 
[3, H, H, H[. Detailed comparisons using direct numeri- 
cal simulations were carried up to large Reynolds num- 
bers and spatial resolutions @ . We can think about this 
problem following an analogy by Betchov 0: a bag full 
of nails has after being shaken its nails pointing in every 
direction in space, and the resulting spatial distribution 
is mirror-symmetric. On the other hand, a well-shaken 
bag full of right-handed screws has its screws pointing in 
any direction, but the mirror-symmetry is broken by the 
screws. In isotropic and homogeneous turbulence, the 
presence of helicity does not changeproperties of turbu- 
lence such as the energy spectrum [J [B|, ||| or the energy 
decay rate @; both the energy and the helicity cascade 
towards smaller scales with the same time scales. 

However, the presence of rotation breaks this degener- 
acy, as the dynamical equations for the flow evolution are 
now sensitive to mirror reflections. Rotation also breaks 
down the isotropy of the flow, introducing a preferred 
direction. The development of anisotropics in a rotating 
flow has been studied in great detail [TO, [0, G3 . How- 
ever, the role played by helicity in these flows has been 
given less attention with only a few exceptions [TH, [TH 
even though the study of helical rotating flows is relevant 
for many atmospheric phenomena [lH llfjl , [l7j | . 

In a recent series of papers [TH, Qj| , we presented evi- 
dence of differences in the scaling laws of rotating flows 
with and without helicity. The presence of helicity in 
a rotating flow changes the energy scaling, as shown in 
numerical simulations and explained with a phenomcno- 
logical theory [lj|| . Changes in the directions of the en- 
ergy and helicity cascades and their associated scaling 
laws have implications for the decay and predictability 



of a helical rotating flow: helical rota ting flows decay 
slower than non-helical rotating flows [20[. These dif- 
ferences imply that to model helical rotating flows in na- 
ture, subgrid models that take into account contributions 
to the turbulent transport coefficients from the helicity 
are required. In agreement with this, a sub-grid scale 
model based on the eddy-damped quasi-normal Marko- 
vian (EDQNM) closure 0, HH, Hl| and that takes into ac- 
count both the cascades of energy and of helicity, proved 
to behave better at reproducing simulations of rotating 
turbulence than models based solely on the energy cas- 
cade im. 

However, the grid resolution of these previous simu- 
lations of rotating helical turbulence (up to 512 3 grid 
points) was insufficient to study together and at suffi- 
ciently high Reynolds number the direct and inverse cas- 
cades; rather, these cascades were studied separately [l9| . 
Also, the amount by which anisotropics develop at small 
scales (or the possible recovery of isotropy at small scale) 
was insufficiently quantified. In order to go further in our 
analysis of rotating turbulence, we present in this paper 
a detailed study of the results of two direct numerical 
simulations of rotating turbulence at unprecedented res- 
olutions. The spatial resolution attained in the simula- 
tions allows us to confirm the scaling laws for the energy 
and helicity spectra predicted in [l9( in runs where the 
direct and inverse cascades now coexist, each with well 
defined inertial ranges. We also analyze the evolution of 
global quantities in the simulations, the development of 
anisotropics in the flow, as well as scaling laws in the 
directions parallel and perpendicular to the axis of rota- 
tion. Section [TT] describes the simulation and the numer- 
ical methods used, Sec. IIIII describes the time evolution 
of global isotropic and anisotropic norms (energy, helic- 
ity, dissipation, characteristic length scales), and Sec. IIVI 
considers the energy and helicity spectra as well as their 
respective fluxes. We also compare with previous runs 
at different Rossby and Reynolds numbers in Sec. [V] 
to identify trends and dependencies with the controlling 
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parameters, the Reynolds and Rossby numbers. Finally, flows [311 ] . 
Sec. IVII gives our conclusions. 

In summary, the results presented in this paper pro- 
vide a detailed description of the numerical simulations 
and compute global and spectral quantities often con- 
sidered in previous studies of rotating turbulence. In a 
following paper (Paper II), we consider in detail the in- 
tcrmittency (or lack thereof) of the energy and helicity 
direct cascades, and study structure functions as well as 
probability density functions of velocity and helicity in- 
crements. Previous studies of rotating turbulence from 
numerical simulations [H, H3 showed a substantial de- 
crease in the intermitt ency of the flow when rotation is 
Hi, IH| for experimental results). 



strong (see also |2£ 
The present simulations can be used to quantify the im- 
pact of helicity in the flow intermittency, and the analysis 
presented in this first paper serves as a reference to allow 
comparisons between these large-resolution simulations 
and previous simulations of rotating turbulence in the 
literature. Overall, these studies allow us to consider the 
recovery of isotropy at small scales, as well as the de- 
velopment of structures at large and small scales in the 
flow, and give a thorough description of rotating turbu- 
lence at scale separations not considered before in direct 
numerical simulations. 



II. NUMERICAL SIMULATIONS 

We solve numerically the equations for an incompress- 
ible rotating fluid, 

— +uj x u + 2f2 x u = -VP + fV 2 u + F, (1) 



and 



k 2 

F = F ^2 {[Bcos(k F y) + Csm(k F z)]x+ 

kp—k\ 

+ [C cos(kpz) + Asm(kpx)} y + 
+ [A cos(kpx) + B sin(kpy)] z} , 



(3) 



where Fq is the forcing amplitude, A = 0.9, B = 1, 
C = 1.1 [32]. An ABC flow, as e.g. in Eq. © for 
only one value of hp, is an eigenfunction of the curl with 
eigenvalue kp; as a result, when used as a forcing func- 
tion, it injects both energy and helicity in the flow. It 
should be noted that in homogeneous turbulence the he- 
licity spectrum cannot develop if it is initially zero (see 
e.g., [Tol. fllj). or if an external mechanism does not in- 
ject helicity. In nature, helicity is created e.g., in the 
presence of rotation and stratification [33|, or near solid 
boundaries in rotating vessels [34| . The use of the ABC 
forcing, although artificial, allows us to study helical ro- 
tating turbulence without the extra computational cost 
associated to the presence of boundaries or stratification. 

Two simulations will be considered in Sees. IIIII and 
HVl with the forcing acting from k\ = 7 to fc2 = 8. This 
leaves some room in spectral space for cascades to develop 
both at large scale and at small scale. The viscosity is 
v = 1.6 x 10~ 4 in both runs, and the time step At = 
2.5 x 10~ 4 . For the first run (run A), ft = 0.06 and 
it is started from a flow initially at rest. The run was 
continued for near 10 turnover times, when a turbulent 
steady state was reached. The value of Fo was such that 

1/2 

in the steady state, the r.m.s. velocity U = (u 2 ) was 
of order unity. For the second run (run B), f2 = 9, and it 
is started from the velocity field in run A at t « 10. Run 
B was continued for 30 turnover times. 

The Reynolds, Rossby, and Ekman numbers of the runs 
quoted in the following sections are defined as usual as: 



V • u = 0, 



(2) 



Re = 



LpU 



(4) 



in a three dimensional box of size 27r with periodic bound- 
ary conditions using a parallel pseudospectral code with 
a spatial resolution of 1536 3 regularly spaced grid points 
(other resolutions will be briefly considered in Sec. fV|) . 
Here u is the velocity field, w = Vxuis the vorticity, and 
v is the kinematic viscosity. The total pressure P mod- 
ified by the centrifugal term is obtained by taking the 
divergence of Eq. ((T|) , using the incomprcssibility condi- 
tion ([2]), and solving the resulting Poisson equation. We 
choose the rotation axis to be in the z direction: fi = Q,z, 
with Q the rotation frequency. Time derivatives are es- 
timated using a second order Runge-Kutta method, and 
the code uses the 2/3-rule for dealiasing. As a result, the 
maximum wave number is fc max = -/V/3 where N is the 
linear resolution. The code is fully parallelized with the 
message passing interface (MPI) library [29|, [3(3 ■ 

The external mechanical forcing F in Eq. |T]) is given 
by a superposition of Arn'old-Bcltrami-Childress (ABC) 



Ro = 



U 



2QL f 



and 



Ek=^ 
Re 



20L 2 ,' 



(5) 



(6) 



where Lp = 27r/min{fci?}, and the turnover time at the 
forcing scale is then defined as T = Lp/U. 

In the following, it will be useful to also introduce a 
micro- Rossby number as the ratio of the r.m.s. vorticity 
to the background vorticity (rotation), 



Ro w = — 



2fi 



(7) 



The value of the micro-Rossby number plays a central 
role in the inhibition of the energy cascade in rotating 
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FIG. 1: Energy (solid), helicity normalized by the forcing 
wavenumber (dashed) , and enstrophy rescaled by the dissipa- 
tion wavenumber ~ 500 (dotted) as a function of time in run 
B with Ro = 0.06. 



FIG. 2: Maximum vorticity in the flow as a function of time, 
for runs A (solid) and B (dotted). 



turbulence 11|. If the micro- Rossby number is too small, 
non-linear interactions are completely damped. Accord- 
ing to (35}, anisotropies develop in rotating flows when 
the Rossby number Ro < I and when the micro-Rossby 
number Ro w > 1 (it is worth noting that the actual values 
for the transition depend on the particular flow studied). 

Based on these definitions, the resulting Reynolds 
number for both runs was Re « 5100. The Rossby num- 
ber of run A is Ro « 8.5, while the Rossby number of 
run B is Ro ps 0.06. This results in Ekman numbers 
Ek w 1.6 x 10~ 3 for run A, and Ek w 1.1 x 10~ 5 for run 
B. The micro-Rossby number of run B is Ro w « 1.2 (in 
all definitions, U and u> are measured in the steady state 
of run A, or when the inverse cascade of energy in run 
B starts). We therefore study flows with large Reynolds 
numbers, but with moderate Rossby numbers as often 
encountered in geophysical problems. Runs A and B will 
also be compared in Sec. fVl with other runs with helical 
forcing at lower resolution as already described in [l9| . 



III. TIME EVOLUTION 
A. Isotropic quantities 

Figure [1] shows the time evolution of the energy, helic- 
ity, and enstrophy in run B. After a transient that lasts a 
few turnover times, energy grows monotonically in time. 
Helicity and enstrophy reach a steady state with nearly 
constant values, smaller than their values at t = which 
correspond to the steady state in run A. 

The enstrophy, proportional to the mean square vortic- 
ity, can be used as a proxy of the relevance of the rotation 
near the Taylor scale, as the ratio of the r.m.s. vorticity 
to the background rotation is proportional to the micro- 
Rossby number. However, turbulence is characterized by 
strong fluctuations of quantities in space and time, and 
one may ask how relevant is the background rotation in 
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FIG. 3: Integral scales associated with the energy (solid line) 
and the helicity (dash line) as a function of time in run B. 



the structures that correspond to these fluctuations. An- 
other measure of the relevance of rotation at the small 
scales can thus be obtained by looking at regions in the 
flow with maximum vorticity. Fig. [2] shows the time his- 
tory of the maximum of vorticity in runs A and B. In 
run A, which starts from a fluid at rest, maxjw} grows 
rapidly from zero and after reaching a maximum sat- 
urates near 800 with strong fluctuations around the 
mean. In run B, started from the last snapshot of run 
A, maxju} starts from the previous value and reaches 
a maximum as the flow becomes anisotropic and the in- 
verse cascade develops, and later saturates with a time 
average value of max{w} 630 with strong peaks. As a 
result, in run B max{uj}/Q « 70. 

The increase of energy with time observed in Fig. [1] 
is associated to an inverse cascade of energy that results 
in an increase in the characteristic size of the energy- 
containing structures in the flow. This is illustrated in 
Fig. [3] by the time evolution of the isotropic integral 
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FIG. 4: Ratio of perpendicular to parallel integral scales, as- 
sociated to the energy (solid) and to the helicity (dashed) as 
a function of time in run B. 
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FIG. 5: Ratio of perpendicular to parallel Taylor scales, as- 
sociated to the energy (solid) and to the helicity (dashed) as 
a function of time in run B at Ro = 0.06. 



scales of the energy 



L 



~E ^ 



E(k) 



k=l 



and of the helicity 



fe=i 



H{k) 
k ' 



(8) 



(9) 



where E and H are respectively the energy and the helic- 
ity, and E(k) and H (k) are the isotropic energy and helic- 
ity spectra. The wavenumber k corresponds to the mean 
radius of the spherical shell containing all modes with 
wave vectors with length between k — 0.5 and fc + 0.5. Ini- 
tially, both integral scales in run B are close to Lp « 0.9. 
However, as time evolves the integral scale of the energy 
increases, while the integral scale of the helicity remains 
approximately constant after a short transient increase. 



B. The development of anisotropy 

The flow in run B becomes anisotropic after a few 
turnover times. Several indicators can be defined to 
quantify the development of anisotropy. In this section, 
we will focus on global quantities, and later we will con- 
sider anisotropy from the point of view of Fourier spectra. 
Directional integral scales as in Eqs. (JSJ) and (0 can be 
defined integrating spectra over wave vectors parallel and 
perpendicular to the rotation axis. We will consider the 
energy spectra E(k\\) and E(k±), where E{k\\) is obtained 
by integrating the energy in all modes with wave vectors 
with k z = fcii between fcii —0.5 and ku +0.5 (i.e., integrat- 
ing over planes in spectral space), and E(k±) is obtained 
by integrating the energy in all modes with wave vectors 
with [k 2 x + fcy) 1 / 2 between kj_ - 0.5 and k± + 0.5 (i.e, 
integrating in cylindrical shells; see e.g., 0, These 
spectra are often referred to in the literature as "reduced" 



energy spectra. Similar procedures can be used to con- 
struct the reduced helicity spectra H(k^) and H(k±). 
Then, following Eqs. j8]) and ((9]), we define perpendicu- 
lar and parallel integral scales 
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(11) 



where the subindcx {_L, ||} denotes that either parallel or 
perpendicular wave vectors are used. 

Figure U shows the perpendicular-to-parallcl ratio of 
integral scales for the energy and the helicity in run B. 
As Q is suddenly increased at t = from its previous 
value in run A, these two ratios first decrease from their 
initial values. Then, / increases slightly and seems 
to reach a steady state after t « 10, while L±/Ln keeps 
increasing monotonically in time on the whole, following 
the increase in the energy due to the inverse cascade. As 
a result, the energy (dominated by large scales) seems 
to become more anisotropic than the helicity which is 
concentrated in smaller scales. 

Integral scales are dominated by the contributions from 
the energy containing scales. Note that both L and L H 
are, after the short transient, larger than the forcing scale 
Lp (sec Fig. [3]). The ratios studied previously then give 
a global indication of anisotropics at scales that are in the 
inverse cascade range. Anisotropies in the (small scale) 
direct cascade range can be quantified, e.g., by the ratio 
of the perpendicular-to-parallel Taylor scales (Fig. O. 
The perpendicular and parallel Taylor scales based on 
the energy and the helicity are defined as 
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FIG. 6: Shebalin angles based on the energy (solid) and he- 
licity spectra (dashed) as a function of time in run B. 



FIG. 7: Ratios E{k\\ = 0)/E (solid line) and H(k\\ = 0)/H 
(dash line) as a function of time in run B. 
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After a transient, the ratio A^/Aj^ in Fig. [5] stabilizes 
at its original value at t = 0, while Aj_/A|| increases with 
time. However, the increase in this ratio is slower than in 
the case of L±/Lu, and presents stronger fluctuations in 
time, being associated with scales in the inertial range of 
the direct energy and helicity cascades. Note that, at the 
onset of the inverse cascade for run B, the perpendicular 
Taylor scale is Aj_ « 0.25, giving for the Taylor Reynolds 
number of that flow (based on the perpendicular scale), 
Ra_l = UX±/v w 1600 (for run A, we have R A « 900 in 
the turbulent steady state; the increase of R> in run B is 
associated to the anisotropization and increase of char- 
acteristic scales of the flow when rotation is increased). 
Note that these values arc larger than the values con- 
sidered in experiments with similar Rossby numbers (see 
e.g., HI). 

Yet another measure of small scale spectral anisotropy 
is given by the Shebalin angles [13] , 



tan 2 (6>) = 2^ 



k 2 ± E(k A 



kf,E(k« 



tan z 



= 2 



£ fcx kjH(k ± ) 



(14) 



(15) 



These angles measure the spectral anisotropy level, with 
the case tan 2 (6) = 2 corresponding to an isotropic flow. 
As the previous quantities, they only give a global mea- 
sure of small-scale anisotropy, and are a byproduct of ax- 
isymmetric energy spectra (see [ll|, H3 ) ■ Figure [5] shows 
the time evolution of the angles based on the energy and 
on the helicity. The helicity at small scales is again more 
isotropic than the energy. However, unlike the previ- 
ous quantities, the Shebalin angles grow fast and then 



saturate in both cases, reaching a steady state after 10 
turnover times. 

Finally, the amount of energy and helicity in two- 
dimensional modes can be measured with the ratios 
E(k\i = 0)/E and H(k\\ = 0)/H (see Fig. EJ). Again, the 
spectral distribution of energy is more anisotropic than 
for helicity. Note that at late times a substantial fraction 
of the energy is in modes with k\\ = 0; at t « 29 near 95% 
of the energy is in those modes, while less than 75% of the 
helicity is in the same modes. All these results indicate 
that the distribution of energy is more anisotropic than 
that of helicity at all scales. As will be discussed next, 
this is due to the fact that helicity only suffers a direct 
cascade and is therefore transported in spectral space to 
smaller scales which are more isotropic. 



IV. SPECTRAL BEHAVIOR 

Figure[8]shows the isotropic energy and helicity spectra 
in run A. The run, with negligible rotation effects, dis- 
plays the usual Kolmogorov scaling in the inertial range 
of the energy and the helicity, with a dual cascade of 
both quantities towards small scales. As in many simu- 
lations of three-dimensional isotropic and homogeneous 
turbulence, the short inertial range is followed by a bot- 
tleneck (which makes the spectra slightly shallower) and 
then by a dissipative range. The dual cascade towards 
smaller scales is further confirmed by examination of the 
energy and helicity fluxes (inset of Fig. |5J) which are 
both positive and constant across the inertial range to 
the right of the forcing wavenumbcr. At wavenumbers 
smaller than fop, both fluxes are negligible. The small 
amount of energy and helicity observed in the spectra at 
those wavenumbers is the result of backscatter, not of 
a cascade, and the energy in the large scales displays a 
slope compatible with a ~ k 2 scaling (see e.g., |38l. |39|). 

The energy and helicity spectra and fluxes at late times 
in run B at Ro = 0.06 arc shown in Fig. [5J An inverse 
cascade of energy develops, as evidenced in the spectrum 
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FIG. 8: Isotropic energy spectrum (solid) and helicity spec- 
trum (dashed) normalized by the forcing wavenumber in run 
A with Ro=8.5. Kolmogorov scaling is shown as a reference. 
The inset gives the isotropic energy flux, and the helicity flux 
normalized by the forcing wavenumber. 



FIG. 10: Perpendicular energy and helicity spectra in run B 
with Ro = 0.06, with same labels as in Fig. [8] Kolmogorov 
and k~ 2 scaling laws are shown as a reference. The inset 
shows the perpendicular energy and helicity fluxes. 
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FIG. 9: Isotropic energy and helicity spectra in run B with 
Ro = 0.06 (same labels as in Fig. [HJ); k~ s/:i (Kolmo gorov), 
and k~ 2 scaling laws are shown as a reference. The inset 
gives the isotropic energy and helicity fluxes. Note the excess 
of (normalized) helicity and of its flux in the small scales. 



by the pile up of energy at scales larger than the forc- 
ing, and in the energy flux by a range of wavenumbcrs 
with nearly constant and negative transfer. However, un- 
like two-dimensional turbulence [ioj . not all the energy 
injected in the system undergoes an inverse cascade: a 
substantial fraction of the injected energy (approximately 
half at this Rossby number) is still transferred to small 
scales in a direct cascade of energy. Moreover, this direct 
cascade of energy is sub-dominant to a direct cascade of 
helicity. All the helicity injected in the system cascades 
to small scales, and as a result the helicity flux (prop- 
erly adimensionalized) is larger than the energy flux at 



FIG. II: Parallel energy and helicity spectra in run B, with 
same labels as in Fig. [8] The inset shows the parallel energy 
and helicity fluxes. 

all wavenumbers larger than Uf- 

Figure [TOl displays the perpendicular energy and helic- 
ity spectra in run B at the same time as in Fig. [9l The 
same features as in the isotropic spectra can be iden- 
tified, and almost no differences arc observed between 
these spectra and the ones showed in Fig. [9l This is 
in agreement with the fact that at late times, ~ 95% 
of the energy and ~ 75% of the helicity are in modes 
with fc|| = (i.e., wave vectors perpendicular to the axis 
of rotation). Here and in the isotropic case, the energy 
spectrum is slightly steeper than the helicity spectrum, 
the result of a dominant direct cascade of helicity. 

No clear scaling is observed in the parallel spectra (Fig. 
[rT]) . There is an excess of helicity at small scales when 
compared with the energy, as in the previous spectra, but 
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FIG. 12: Product of the energy and helicity spectra in run 
B, normalized by &f and compensated by k 4 . The solid line 
corresponds to the isotropic spectra, the dash line to the per- 
pendicular spectra, and the dotted line to the parallel spectra. 



FIG. 13: Ratio (AE + AI1)/AE as a function of the in- 
verse Rossby number for several runs (see text for definition). 
Squares correspond to runs A and B. Diamonds correspond 
to helical runs with 512 3 resolution, Uf = 7, and Re « 1200 
for different Rossby numbers, as analyzed in 



here an inertial range cannot be identified. Moreover, the 
energy and helicity fluxes in the parallel wave vectors are 
positive at all scales, and show no approximately con- 
stant range. The fluxes peak at kn = 1 and decrease fast 
for larger wavenumbcrs. This indicates that in the par- 
allel direction, a small portion of the energy transferred 
towards large scales by the inverse cascade of energy in 

, is transferred back towards small scales although not 
through a cascade. The wavenumber band where forcing 
occurs is not identifiable either. 

The presence of waves in a rotating flow slows down 
the direct energy cascade, resulting in a steeper energy 
spectrum than in isotropic and homogeneous turbulence. 
However, unlike isotropic and homogeneous turbulence, 
the slope of the energy spectrum depends on whether the 
flow is helical or not. This is the result of the flow de- 
veloping, at small scales, a dominant direct cascade of 
energy in the former case, and a dominant direct cas- 
cade of helicity in the latter. In a flow with maximum 
helicity (H(k) = kE(k)), the energy and helicity spectra 
are, according to these arguments (see [lj| for a detailed 
derivation), predicted to be E ~ kj 2 ' 5 and H ~ fcj 1 " 5 . 
Note that this results in the relative helicity H(k)/[kE(k) 
independent of wavenumber, i.e. corresponding to an 
alignment of velocity and vorticity identical at all scales. 
In the general case (not maximally helical), the energy 
spectrum gets closer to a kj 2 spectral law, but with both 
spectra such that their product is still E(k)H(k) ~ fcj 4 . 
The slope of the energy spectrum in the direct inertial 
range of Figs. 191 and [TfJl is « 2.1, while the slope of the 
helicity spectrum is w 1.9, in good agreement with these 
predictions (slopes where obtained through a least square 
fit). Figure [12] shows the product of all spectra compen- 
sated by k . The isotropic and perpendicular spectra 
show a flat region compatible with this scaling, and the 
results are in good agreement with the arguments based 
on a dominant cascade of helicity in the helical rotating 
case. No clear scaling is observed in the parallel direction. 



V. SCALING WITH REYNOLDS AND ROSSBY 
NUMBERS 

In this section we analyze an ensemble of runs in or- 
der to see the emergence of scaling laws when examin- 
ing a range of Reynolds and Rossby numbers; to that 
effect, we combine runs A and B of this paper with pre- 
vious simulations at different Reynolds and Rossby num- 
bers fl9l ] . The sub-dominant direct cascade of energy in 
run B is the result of the relation e ~ kpe at the forc- 
ing scale (where e is the helicity injection rate dH/dt 
and e the energy injection rate dE/dt), together with 
the development of an inverse cascade of energy at small 
enough Rossby number which removes some of the in- 
jected energy that would otherwise be available for the 
direct cascade. Since a fraction of the energy injected at 
kp goes towards large scales (in run B, | H | w e/2, where 
II - is the energy flux towards large scales, measured at 
wavenumbers smaller than kp), while most of the helicity 
injected goes towards small scales (with E + the helicity 
flux towards small scales) , the energy flux towards small 
scales n + is smaller than E + in the entire direct cas- 
cade range (see e.g., Figs. l9l and [T0|) . This results in a 
direct cascade dominated by the helicity, reminiscent of 
the pure direct cascade of helicity hypothesized in [4lj . 
However, here the time scale of the direct cascade is af- 
fected by the presence of Rossby waves, resulting in the 
E(k)H(k) ~ fcj 4 rule as discussed in [lj|. 

This scenario is mainly based on two hypothesis, that 
the injection rates are related through e ~ kpe (with the 
proportionality constant equal to one when the forcing 
injects maximum helicity at only one wavenumber, and 
smaller than one when the forcing is not maximally he- 
lical), and that the ratio E + /n + in the direct cascade 
range increases with rotation (cither monotonically or 
saturating at a small value of Ro such that the direct 
helicity flux becomes dominant). 
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FIG. 14: Ratio £ + /(fcFiI + ) as a function of the inverse 
Rossby number for several runs (see text for definitions). 
Squares correspond to runs A and B. Diamonds correspond to 
helical runs with 512 3 resolution, Uf = 7, and Re « 1200. Tri- 
angles correspond to helical runs with 512 3 resolution, fci? = 2, 
and Re « 5700 0. 

That the former condition holds, in the range of pa- 
rameters tested in the present study, is illustrated in Fig. 
IT3"1 which displays the normalized ratio (AE+AII)/A£ as 
a function of the Rossby number for several runs, where 

AS = (£+ - TT)/k F = e/k F , 

and 

An = n+ - n - = e . 

Note that these differences (the amount of flux to small 
scales plus the amount of flux to large scales, since the 
latter is negative) are proportional to the helicity and 
energy injection rates e and e. For a forcing that injects 
maximum helicity, if e = kpe, then (AE + AF1)/AE = 2. 
The constancy of (AE + AII)/AE for runs at different 
Reynolds and Rossby numbers, as observed in Fig. 1131 
builds confidence on the validity of the former hypothesis. 

Figure [Til shows the ratio E + /(fci?II + ), i.e., the ratio 
of energy flux towards small scales to the normalized he- 
licity flux towards small scales, as a function of Ro -1 . 
While for large values of the Rossby number this ratio is 
close to unity (where a dual cascade of energy and helicity 
takes place, with both quantities having the same spec- 
tral index in the inertial range), as the Rossby number 
is decreased the ratio becomes larger than one. Simu- 
lations with very different Reynolds numbers and scale 
separation between the domain size and the forcing scale 
appear to collapse onto one curve that grows as Ro de- 
creases, thus building confidence on the validity of the 
latter hypothesis. 

VI. CONCLUSIONS 

We presented results from two simulations of rotat- 
ing turbulence at large Reynolds numbers and moderate 



Rossby numbers, with spatial resolution of 1536 3 grid 
points and forced at intermediate scales, allowing for si- 
multaneous direct and inverse cascades of the ideal in- 
variants of the flow to develop. The forcing injects both 
energy and helicity, and both inverse and direct cascades 
of energy develop in the case of stronger rotation, to- 
gether with a direct cascade of helicity. The inverse cas- 
cade range is dominated by the energy, and the flow is 
anisotropic and only weakly helical. On the other hand, 
in the direct cascade range, the normalized helicity flux 
is larger than the energy flux. The dominance of the he- 
licity flux in this range results in the time scales of the 
direct cascade being imposed by the helicity. As a result, 
the small scale helicity spectrum is shallower than the 
energy spectrum. Both spectra scale as E rj kT^ and 
H « fcj with e + h = 4 as predicted by phenomenolog- 
ical arguments (for the run with Ro » 0.06, it was 
found e w 2.1 and h « 1.9). The slope of the energy 
spectrum is slightly steeper than what is found in recent 
numerical simulations of non-helical rotating turbulence 
(see e.g. MM)- 

Comparisons with other simulations of helical rotat- 
ing flows albeit at lower resolution allowed us to build 
confidence on the dominance of the direct cascade of he- 
licity over the energy as the Rossby number is decreased. 
While direct and inverse inertial ranges of energy and 
helicity were identified in the isotropic and perpendicu- 
lar spectra and fluxes, no clear scaling was found in the 
parallel direction. The energy and helicity fluxes in this 
direction are positive for all wavenumbers, and no range 
of scales with constant flux was found. 

The development of anisotropies in the flow was stud- 
ied using global and spectral quantities. In all cases, 
it was observed that the distribution of helicity is more 
isotropic than the distribution of energy. As an example, 
it was found that at late times w 95% of the energy is in 
modes with kn = 0, while less than 75% of the helicity is 
in the same modes (as a comparision, for run A « 24% 
of both the total energy and helicity are in modes with 
fen = 0). Other measures of anisotropy at both large and 
small scales gave consistent results. The vortical struc- 
tures that develop in such flows and the anisotropy, as 
well as the recovery of isotropy at small scales, will be 
studied in more detail in Paper II using probability den- 
sity functions of velocity and helicity increments, struc- 
ture functions based on the symmetries of the problem 
and three-dimensional visualizations of the velocity and 
vorticity. 
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